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With Hubbard model, the entanglement scaling behavior in a two-dimensional itinerant system 
is investigated. It has been found that, on the two sides of the critical point denoting an inherent 
quantum phase transition (QPT), the entanglement follows different scalings with the size just as 
an order parameter does. This fact reveals the subtle role played by the entanglement in QPT as a 
fungible physical resource. 
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o 

The existence of entanglement between distinctive quantum systems has marked a fundamental difference between 
quantum and classical physics. Recently, with the explosive development of the research in quantum information 
£>V theory and quantum computation 0, 0, 3 , the study of the entanglement 0, 0] has come to the limelight again 
after more than 60 years of controversies and strenuous progress. Experimentally, entanglements have already been 
produced between up to four photons 0, @ and even between two macroscopic states such as two superconducting 
qubits, each of which contains as large as 10 9 electrons 0. But theoretically, because an ensemble's Hilbert space 
grows exponentially with the number of its component particles, we are still far from fully understanding the contents 
of the entanglements. Only for the simplest state with two distinguishable particles can we have a complete description 
I . of the entanglement measure. For states of more than two particles, especially for mixed states, the current knowledge 
' about their entanglement is very limited and all the related complexities have just begun to be explored. For spin-only 
[*""». , entanglement of localized distinguishable particles, the most popular measure of the entanglement is the Wootters' 
measure 10]. Recently, the influence of the quantum statistics on the definition of the entanglement has begun to 
be noticed and discussed by several authors |lll Il2l . Although various entanglement measures have been put 
forward, according to Gitting's criterions ^5|, only Zanardi's measure ^(| survives the test of all the requirements 
upon entanglement definition. This measure is given in Fock space as the von Neuman entropy, namely, 

o : 

^ ' Ej = -Trpj Inpj, pj = Trj \tp) (tp\ , (1) 

where Trj denotes the trace over all but the jth site and ip is the antisymmetric wave function of the studied system. 
Hence Ej actually describes the entanglement of the jth site with the remaining sites. A generalization of this one-site 
entanglement is to define an entanglement between one L-site block with the rest of the systems ^3] , 

cr 



E L = -Tr{p L \og 2 p L ). (2) 



where all the sites are traced out except those belonging to the selected block. 

Recently, it has been speculated that the most entangled systems could be found at the critical point when the 
system undergoes a quantum phase transition (QPT), i. e. a qualitative change of some physical properties takes 
place as an order parameter in the Hamiltonian is tuned |l8j . QPT results from quantum fluctuations at the absolute 
zero of temperature and is a pure quantum effect featured by long-range correlations. So far, there have already 
been some efforts in exploring the above speculations, such as the analysis of the XY model about the single-spin 
entropies and two-spin quantum correlations [l9l I20I] , the entanglement between a block of L contiguous sites and 
the rest of the chain and also the scaling of entanglement near QPT [2l| . But because there is still no analytical 
proof, the role played by the entanglement in quantum critical phenomena remains elusive. Generally speaking, there 
exists at least two difficulties in resolving this issue. First, until now, only two-particle entanglement is well explored. 
How to quantify the multi-particle entanglement is not clear. Second, QPT closely relates to the notorious many- 
body problems, which is almost intractable analytically. Until now, the only effective and accurate way to deal with 
QPT in critical region is the density-matrix renormalization group method |22|. Unfortunately, it is only efficient for 
one-dimensional cases because of the much more complicated boundary conditions for two-dimensional situation. It 
should be mentioned here that recently, Vial [231 ] has put forward another new efficient numerical method to study 
one-dimensional many-body systems based upon the entanglement contained in the system. 

In this paper, we will focus on investigating the entanglement behavior in QPT for two-dimensional array of 
quantum dots, which provide a suitable arena for implementation of quantum computation |24l l25l 1,26]. For this 
purpose, the real-space renormalization group technique |27| will be utilized and developed for the finite-size analysis 
of entanglement. 
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The model we use is the Hubbard model with the Hamiltonian, 

H= -t J2 [4,c j<r + H.c] 

<i,j>,a 

+UY,i\-^){\-n il )+KY,Ii (3) 

i i 

where t is the nearest-neighbor hopping term, U is the local repulsive interaction and K — — U/4 and Ii is the unit 
operator. c^.(cj CT ) creates (annihilates) an electron with spin a in a Wannier orbital located at site i; the corresponding 
number operator is rii a = c^.Cj CT and <> denotes the nearest-neighbor pairs. H.c. denotes the Hermitian conjugate. 

For a half-filled triangular quantum lattice, there exists a metal-insulator phase transition with the tuning parameter 
U/t at the critical point 12.5 psl 12^. l30| . The corresponding order parameter for metal-insulator transition is the 
charge gap defined by A g = E(N e — l) + E(N e + l) — 2E(N e ), where E(N e ) denotes the lowest energy for a N e — electron 
system. In our case, N e is equal to the site number N s of the lattice. Unlike the charge gap calculated from the 
energy levels, the Zanardi's measure of the entanglement is defined upon the wave function corresponding to E(N e ) 
instead. Using the conventional renormalization group method for the finite-size scaling analysis [28l I2II l30l| . we can 
discuss three schemes of entanglement scaling: 

1) Single-site entanglement scaling with the total system size, E S i ng i e ; 

2) Single-block entanglement scaling with the block size, Euock] and 

3) Block-block entanglement scaling with the block size, Euock-biock- 

Fig. 1 presents the single-site entanglement scaling. It is obvious that E s i ng i e is not a universal quantity. This 
conclusion is well consistent with the argument given by Osborne |19|, who claims that the single-site entanglement 
is not scalable because it does not own the proper extensivity and does not distinguish the local and distributed 
entanglement. One more interesting feature in Fig. 1 is that when the system size is increased beyond 7 2 , E S i ng i e 
almost makes no change any more. 

This implies that only a limited region of sites around the central site contributed significantly to the single-site 
entanglement. Using the one-parameter scaling theory, near the phase transition point, we assume the existence of 
scaling function / for Euock-Uock such that: 

Euock-biock =q VE f(j) (4) 

where q = U/t—(U/t) c measures the deviation distance of the system away from the critical state with (U/t) c — 12.5, 
which is exactly equal to the critical value for metal- insulator transition when the same order parameter U/t is used 
[lil |2^, ■ £ = q~ v is the correlation length of the system with the critical exponent v. 
Hence, 

E b iock-block = q yE f{N^q), (5) 

where we used N = L 2 for the two-dimensional systems. 

In Fig. 2, we show the results of Euock-Uock as a function of (U /t) for different system sizes. With proper scaling, 
all the curves collapse onto one curve, which can be expressed as 

Euock-biock = f(qN*). (6) 

Thus the critical exponents in Eq. (5) are Ue = 0, v = 1. It is interesting to note that we obtained the same v as in 
the study of metal-insulator transition. This shows the consistency of the results since the critical exponent v is only 
dependent on the inherent symmetry and dimension of the investigated system. 

Another significant result lies in the finding that the metal state is highly entangled while the insulating state, is 
only partly entangled. For a 4 dimensional density matrix, the maximally entangled state can be written as a diagonal 

4 

matrix with equal components \. The related entanglement is — ^2 1 1°S2 \ = 2, which is exactly the value obtained 

i=l 

from Fig. 2(b). However, unlike the metal state, the insulating states should be expected to have electrons showing 

less mobility. If we assume the highly probable situation, i.e. the central site has equal probability to be in ||), \l) and 

3 

no occupation in |0), |t|), the corresponding entanglement is then Ebiock-block — ~ 2 \ l°g2 4 = 1) a ^ so consistent 

»=2 

with the results from Fig. 2(b). 

All the above discussions are confined to the entanglement between the central block and its surrounding blocks. 
Because the central block is a very special one showing the highest symmetry, one may wonder what can happen to the 
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neighboring blocks, for example, the entanglement between block 7 and the rest 6 blocks. To answer this question, the 
same calculations are conducted and the results are the same except that in the metal state, the maximal entanglement 
is a little less than 2 and the minimal one is a little less than 1. This can be explained by the asymmetric position of 
site 1 in the block. 

It should be mentioned that the calculated entanglement here has a corresponding critical exponent he = 0. This 
means that the entanglement is constant at the critical point over all sizes of the system. But it is not a constant 
over all values of U ft. There is an abrupt jump across the critical point as L — > oo. 

If we divide the regime of the order parameter into non-critical regime and critical regime, the results can be 
summarized as follows: In the non-critical regime, i.e. U/t is away from (U/t) c , as L increases, the entanglement will 
saturate onto two different values depending on the sign of U/t — (U/t) c ; At the critical point, the entanglement is 
actually a constant independent of the size L. 

These properties are qualitatively different from the single-site entanglement discussed by Osborne [l9j , where the 
entanglement with Zanardi's measure increases from zero to the maximum at the critical point and then decreases 
again to zero as the order parameter 7 for XY mode is tuned. 

These peculiar properties of the entanglement we have found here can be of potential interest to make an effective 
ideal "entanglement switch". For example, with seven blocks of quantum dots on triangular lattice, the entanglement 
among the blocks can be regulated as "0" or "1" almost immediately once the tuning parameter U/t crosses the 
critical point. The switch errors will depend on the size of the blocks. Since it has already been a well-developed 
technique to change U/t for quantum dot lattice |2y,|3j, the above scheme should be workable. To remove the special 
confinement we have made upon the calculated entanglement, namely only the entanglement of block 1 and block 7 
with the rest ones are considered, in the following, we will prove that the average pairwise entanglement also has the 
properties shown in Fig. 2(b). 

Let E{,iock~biock,7 denotes the entanglement between the 7th and all the remaining blocks in a hexagonal sys- 
tem. From the symmetry of the system we can show that the total pairwise entanglement E to t — QEuock-Uock + 
3Euock-biock,7- The average 2-site entanglement is E average = E tot /21 = (2Ebiock-Uock + Eu ck-biock,i) /7 '. Because 

E b iock,7 = 9i(qN*), Euock,i = 9a(qNi), (7) 

then we should have E average = g(qN 1 ^ 2 ), where 51,(72 and g are scaling functions. 

For obtaining the single-block entanglement, the first step is to make a cutoff over the system size. In our work, we 
let it be 7 9 . The results are presented in Fig. 2(c). It is magnificent that as we change the size of the central block, 
its entanglement with all the rest sites follows the same scaling properties as Euock-block- It is understandable if we 
consider the fact that only a limited region round the block contribute mostly to Ebiock- This result greatly facilitate 
the fabrication of realistic entanglement control devices, such as quantum gates for quantum computer, since we don't 
need to delicately care about the number of component blocks in fear that the next neighboring or the next-next 
neighboring quantum dots should influence the switching effect. 

In summary, in this paper, various schemes of the finite-size scaling properties with Zanardi's measure of entangle- 
ment has been investigated for Hubbard model on a triangular quantum lattice. The critical exponent v — 1 has been 
found, which coincides well with our previous work in studying a quite different physical property, the charge gap. 
When the block size L — > 00, the entanglement shows an abrupt change when the tuning parameter crosses the phase 
transition point. This property might be well applied to make an "entanglement switch" and shows the promising 
prospect of regarding entanglement as a new physical resource. 
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FIG. 1: Figure 1: (a) Schematic diagram showing the central site (purple) and the surrounding ones in the triangular quantum 
dot lattice. The dotted lines represent the site-site interactions, (b) Scaling of the single-site entanglement for various system 
size. The sizes are denoted by different symbols. 
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FIG. 2: Figure 2 : (a) Schematic diagram displays the lattice configuration with central purple block and the surrounding green 
ones, (b) Scaling of block-block for various system size and (c) Scaling of block entanglements with the block size. 



